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APPLIED MATHEMATICS 


f eee National Bureau of Standards conducts a 
broad program in applied mathematics, oriented 
primarily toward the physical sciences and engineering 
to meet its own needs, but serving also to provide 
mathematical assistance to other Government agencies 
in these and other fields. Although the areas of 
numerical analysis, numerical computation, statistical 
engineering, and mathematical physics receive primary 
emphasis, the range of mathematical applications is as 
extensive as the varied interests of the Bureau’s staff 
and the Government agencies requesting assistance. 
For example, mathematical services are rendered in 
such diverse fields as biology, operations research, 
language translation by electronic computer, and 
management problems. This special issue of the 
Technical News Bulletin indicates the scope of the 
Bureau’s activity in applied mathematics by reporting 
a representative sampling of current work in this area. 

In the laboratories of the Bureau and other Govern- 
ment agencies, modern techniques of statistical in- 
ference play an important role, both in the design and 
the interpretation of experiments. Individual re- 
search workers are becoming increasingly aware that 
consultation with statisticians enables them to decrease 
the cost of experimentation, and, at the same time, to 
Gnerease the amount of information obtained. More- 
over, the statistician’s assistance is essential in the 
planning and evaluation of experimental programs in- 
volving the cooperation of many laboratories, like the 
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interlaboratory standards comparisons in which the 
Bureau participates. Statistical engineering thus per- 
vades the Bureau’s entire program and is a major sub- 
ject of advisory services to other agencies. Typical 
statistical services in the design of experiments and 
the analysis of the reliability of electronic components 
are discussed on pages 20 and 24 respectively. 

Research in numerical analysis is aimed primarily 
at devising computational techniques adapted to high- 
speed electronic digital computers. Sutcessive ap- 
proximation procedures, rather than the use of closed- 
form analytical expressions and sets of tables, are re- 
quired by the nature of these machines. As a numer- 
ical result may stem from thousands of intermediate 
steps, increased attention must be given to the propa- 
gation and accumulation of errors to insure sufficient 
accuracy in the final solution. Also, in numerical 
analysis, it is necessary to devise new  successive- 
approximation techniques suitable for electronic com- 
puters, to investigate their rates of convergence, and to 
determine bounds to the error which may occur in the 
final result. Supporting research is required in a 
broad variety of related fields, such as matrix theory, 
combinatorial analysis, linear programing, game 
theory, and conformal mapping. Some of the work 
in numerical analysis is discussed on page 18. 

The Bureau was one of the pioneers in the art of 
high-speed numerical computation, having designed 
and placed into operation two electronic computers, 
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the SEAC*? and the SWAC,? long before either the 
tremendous potential of such equipment or the feasi- 
bility of its construction was widely recognized. Both 
machines are still productive but the Bureau now op- 
erates extremely fast commercial computational facil- 
ities also. Typical problems programed for its own 
staff and other Government agencies are discussed 
under the headings “Automatic Computation as an 
Aid to Scientific Research”, “High-Speed Computers 
in Radio Propagation Research”, and “Data Process- 
ing and Operations Analysis on a Scientific Com- 
puter’, beginning on pages 2, 10, and 14, respectively. 

The Bureau has collaborated with other groups, both 
nationally and internationally, in the development of 
improved techniques of preparing problems for ma- 
chine solution. The ideal method of problem prepa- 
ration, in both the scientific and business fields, would 
allow the problem to be described by its originator 
in the terminology of his specialty. Although this 
ideal seems far out of reach, at least within the fore- 
seeable future, important steps have been taken in 
devising languages for problem description which are 
easier for the programer to use. Such advances have 


A heat! use of automatic computation as an aid in re- 
search is more and more becoming the rule rather 
than the exception, a requirement rather than an 
option. This growing indispensability of the high- 
speed computer to scientific and engineering research 
is clearly shown by the variety, extent, and importance 
of the problems that are being solved by the use of 
automatic computation at the Bureau. Large-scale 
digital computers help the Bureau to carry out its 
fundamental program of research on properties of 
matter and methods of measurement with increasing 
speed, efficiency, and thoroughness. In addition, they 
aid the Bureau in meeting the scientific computing re- 
quirements of the many Government agencies which 
rely upon NBS experience and facilities in this area. 
The solving of problems arising in science usually 
involves several steps. First, a theoretical analysis is 
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been made possible by experienced programers who 
have devised techniques that direct the computer to 
translate the user-oriented language into appropriate 
sequences of its built-in operations. Work of this 
kind is discussed under “Automatic Programing 
Studies” on page 6, and “OMNIFORM” on page 8. 
The research in applied mathematics also includes 
basic studies conducted to provide mathematical the- 
ories for currently developing scientific and engineer- 
ing fields. This work includes the derivation of solu- 
tions to specific problems either in analytic form or 
in a form suitable for automatic computation. — Illus- 
trations, discussed on pages 25 and 26, are research 
programs in plasma dynamics and satellite trajectory 
theory. Also, extensive fundamental research in linear 
programing and game theory underlies the machine 
simulation studies described on pages 14 and 19. 


*SEAC—the National Bureau of Standards Eastern 
Automatic Computer, NBS Tech. News Bul. 34, 121 
(1950) . 

* National Bureau of Standards Western Automatic 
Computer, NBS Tech. News Bul. 37, 145 (1953). SWAC 
was transferred to the University of California at Los 
Angeles in 1954. 
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Top, left: Artist’s conception of the interaction between 
an atom (white ball) and a diatomic molecule. The 
atom, upon contact with the molecule, loses some of its 
kinetic energy. This energy is transferred to vibrational 
energy, raising the molecule to an excited vibrational 
state. The probabilities for various states as a function 
of the atom’s kinetic energy, as calculated on a high- 
speed computer, are shown above. The 0 > 0 curve rep- 
resents the probability of no vibration, the 0 > 1 curve 
the probability of the molecule being raised to the lowest 
excited vibrational state, and so on. 
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necessary to determine explicit solutions, if they exist. 
Then the best method for numerical solution must be 
selected and, where machine methods are found de- 
sirable, the problem must be programed and coded. 
Sometimes this is all that is required—for example, in 
the setting up of problems for machines of other Gov- 
ernment agencies. This type of service is frequently 
rendered because of the Bureau’s familiarity with 
machine methods of computation. However, for its 
own laboratories and for outside agencies depending 
on the Bureau’s computing facilities, nonroutine prob- 
lems are usually carried to completion on a high-speed 
digital computer, which is supplemented by punched 
card machines and desk computers. 

Work for other agencies constitutes about three- 
fourths of the Bureau’s activities in scientific compu- 
tation. An example of this outside work is assistance 
to the Weather Bureau. The cooperation between 
these two agencies is so extensive that a microwave 
link has been installed between the Weather Bureau’s 
Suitland (Md.) installation and the NBS Computation 
Laboratory. This link, which transcribes data from 
a magnetic tape at either installation to a tape at the 
other, will be especially useful should the Weather 
Bureau’s computation facility ever break down. The 
NBS machine can serve as a standby, making it pos- 
sible to put out weather predictions without interrup- 
tion for time-consuming repairs. Furthermore, the 
link can also save time and effort in transmitting prob- 
lems which the Weather Bureau frequently asks the 
Bureau of Standards to solve. These less-routine com- 
putations demand the greater memory capacity that is 
available at the Bureau’s installation. One such prob- 
lem was an exploratory long-range weather forecast 
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Three-dimensional plot of voltage, distance, and current 
density for argon at a pressure of 10 mm. This model 
is based on vyoltage-current static characteristics calcu- 
lated on a high-speed computer. The work was carried 
out for the Diamond Ordnance Fuze Laboratories, to pro- 
vide an improved basis for gas tube design. Photo 
courtesy of DOFL. 
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made last spring. The machine was used to predict 
atmospheric pressure distribution as long as 2 months 
in advance. 

More typical of services to other agencies is some 
recent work for the National Aeronautics and Space 
Administration and for the Navy’s Vanguard project. 
This included the large-scale reduction of satellite ob- 
servations and predictions of satellite orbits. Other 
NASA problems handled at the Bureau are the com- 
putation of the electrical charge density surrounding 
a satellite, and the solutions of a heat transfer equation 
representing the effects of the Van Allen radiation belt 
over the polar regions of the earth. These problems, 
along with those outlined below, illustrate not only the 
increasing number of fields in which the potentialities 
of the digital computer are being exploited, but also 
the wide range of computational techniques now 
available. 


Gas Tube Characteristics 


In cooperation with the Diamond Ordnance Fuze 
Laboratories, Department of the Army, the Bureau is 
studying the effect of space charge in cold-cathode gas 
tubes. This work involves computations of the volt- 
age-current characteristics of the tubes from the basic 
properties of the gas and electrode material and a 
knowledge of tube geometry. It is hoped that these 
studies will provide a basis for designing tubes for 
optimum performance, eliminating the endless experi- 
mental testing now required. 

The computations, formulated by A. L. Ward of 
DOFL, were programed by I. Stegun, W. Cahill, and 
C. Wade of the NBS Computation Laboratory, in an 
effort to develop a theory for the critical negative-slope 
region which determined the tube’s breakdown char- 
acteristics. Although good agreement previously 
existed between theory and experiment for the very 
low current, positive-slope region and_ satisfactory 
agreement was obtained in the glow region after 
breakdown, no successful theory was available for the 
negative-slope region. 

Preliminary calculations on SEAC (Standards Elec- 
tronic Automatic Computer) at the Bureau showed that 
the effect of space charge upon the electric field in the 
gas can explain this region. Computations were then 
carried out which agreed with experimental measure- 
ments indicating that the space-charge field is a linear 
function of the distance from the cathode to the anode. 
Later computations included the effect of electron space 
charge; the charge density distribution then indicated 
the presence of a positive column plasma. The voltage- 
current static characteristics of all rare gases and of 
some molecular gases have been calculated for parallel 
plate geometries for current values well up into the 
glow discharge region. 


P-Wave Equations 


In working out a theoretical formulation for various 
atomic and nuclear processes, the radiation theory 
eroup has encountered problems which demand 
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machine treatment. A recent problem of this type 
arose in connection with the scattering of electrons by 
hydrogen atoms. This study conducted by A. Temkin 
involved a determination of a phase shift which is 
related to the solution of a differential equation. 
Depending upon the desired degree of approximation in 
the description of the physical process, various differ- 
ential or integro-differential equations can be used. 
They are linear and of the second order with variable 
coefficients which, in the more realistic cases, involve 
integrals of the unknown wave function itself. These 
integro-differential equations can be reduced to a set 
of uncoupled differential equations which, however, 
still require machine computation. A program has 
been written by J. Lamkin, of the radiation theory 
group and about half the desired results have already 
been obtained. Computation of each phase shift re- 
quires less than 1 minute of machine time. 


Eigenvalues of Differential Operators 


An interesting and important property of certain 
linear differential equations with a parameter A is that 
nonzero solutions exist for only a limited number of 
values of A. These values of X are known as the 
“eigenvalues” of the given differential operator. Ina 
physical system, such as a vibrating plate or a rotating 
shaft, they are the frequencies of the various overtones 
of the plate or the critical angular velocities at which 
the shaft will buckle. In the case of the Schroedinger 
equation for an atom or molecule, they represent the 
permissible energies of the system. 

The experimentally observed numbers for such 
eigenvalues often compare remarkably well with the 
numbers calculated as the solution of the theoretically 
derived equations. However, in the majority of prob- 
lems no exact solutions can be found for the eigen- 
values. Therefore, the only approach is to devise tech- 
niques to find numerical estimates. The main idea of 
these approximation methods is to reduce the given 
differential problem to an approximating algebraic 
problem. This is partially accomplished by the classi- 
cal Rayleigh-Ritz procedure, which gives upper bounds 
to the eigenvalues of the differential operator as the 
roots of a secular equation. 


OCHS Deslenations 
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A method for finding lower bounds to the eigen- 
values of differential operators has recently been de- 
vised by N. W. Bazley of the numerical analysis staff. 
This method is a modification of an earlier procedure 
developed by A. Weinstein and N. Aronszajn. It 
applies to cases where the given differential operator 
can be adjusted so that a solvable problem—one which 
gives rough lower bounds to the eigenvalues of the 
original problem—results. This solvable problem is 
linked to the given problem by a sequence of “inter- 
mediate problems” which improve the lower bounds. 
The solution of the intermediate problem reduces the 
determination of the roots of a secular equation. By 
combining the lower bounds with the Rayleigh-Ritz 
upper bounds, one can confine each eigenvalue between 
upper and lower bounds whose difference can be made 
arbitrarily small. The method is being applied to 
several problems of mathematical physics, especially 
the wave equation for helium. The necessary calcula- 
tions are being programed for automatic computation 


Microwave Absorption in Ferrites 


In other work being done for the Diamond Ord- 
nance Fuze Laboratories, the Bureau is evaluating a 
mathematical formulation which predicts the micro- 
wave absorption in ferrites. ‘These magnetic materials 
are used in microwave devices such as isolators, rota- 
tors, and switches. In efforts to extend their use, a 
theoretical prediction, which takes into account re- 
cently observed complicated resonance curves, has 
been derived by O. Cruzan, of DOFL. Such a pre- 
diction, once its results are confirmed experimentally, 
can be used to determine the behavior and properties 
of ferrites with a limited amount of testing, thus aiding 
in the effective utilization of existing materials and 
the synthesis of new materials. 

However, the working out of the proposed formula- 
tion is a tedious and lengthy operation when performed 
by ordinary computation methods. Complex Legendre 
functions are involved. As there are no tables of these 
functions, a complete evaluation must be made. The 
great advantage a computer offers in this connection is 
speed. Evaluating only five points for a function took 
a mathematician working by hand 5 months; with a 


In studying color fading, automatic 
computer techniques are used to 
convert spectrophotometric data into 
Munsell coordinates. J. Schleter 
(left) and H. J. Keegan (standing) 
of the photometry and colorimetry 
group discuss computations with W. 


Rheinboldt (center) and J.. P. 
Menard (right) of the Computation 
Laboratory. 
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machine code written by R. Zucker of the Bureau’s 
Computation Laboratory, 12 points were evaluated in 
31 minutes. 


Analysis of Color Fading 


The program in color analysis has been strongly 
supported by computer activities during the past few 
years. Recently emphasis has shifted to the problem 
of color fading. A determination of fading involves 
first of all the conversion of a large amount of available 
spectrophotometric data from one system to another 
and secondly the computation of the color differences. 
A computer program designed by W. Rheinboldt and 
J. P. Menard of the Computation Laboratory for H. J. 
Keegan of the photometry and colorimetry group 
makes these calculations economically feasible. 

The Bureau’s work requires that the standard CIE- 
chromaticity coordinates for describing color spectro- 
photometrically be translated into the more widely 
used Munsell system, where color coordinates are 
based on perceived value, hue, and chroma. As no 
usable mathematical formula is known for the relation 
between the two systems, the conversion must be based 
on direct experimental comparison of spectrophotom- 
eter measurements with estimates of color spacing by 
human observers. For purposes of computer calcula- 
tion, the results of such comparisons (by a subcom- 
mittee of the Optical Society of America) for about 
5,000 discrete points which form a grid in color space 
are being punched on cards. For colors correspond- 
ing to points other than those on the grid, the con- 
version will be made by a three-dimensional inter- 
polation. 

Converting CIE-chromaticity coordinates for a 
single color to Munsell coordinates and checking the 
results requires about 30 min by hand calculation. A 
computer program, written in 1957, performed the 
same operation in 24 sec. However, this program was 
still too slow for present needs. It slows down the 
machine by requiring, among other things, that all the 
data be stored on magnetic tape. 

More recently a revised program was completed 
which uses the expanded internal memory of the 
machine instead of magnetic tape for storing data. A 
single conversion is now completed in 314 secs. The 
new program also computes color differences from 
Godlove’s formula, where the parameters are Munsell 
renotations of the given samples. During the first 
year of the operation of this program, about 4,600 color 
samples were evaluated. 


Dielectric Relaxation in Gases 


Work on dielectric relaxation of molecules in gases 
is being carried out by H. Oser of the Computation 
Laboratory for the chemistry laboratories. The pur- 
pose of this study is to evaluate resonant and non- 
resonant microwave line shapes of symmetrical-top and 
linear molecules for certain models of relaxation and 
a reasonable range of relaxation times. By comparing 
the results obtained for the proposed models with 
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S. J. Kryder of the chemistry laboratories tunes a resonant 
cavity containing a gas sample in obiaining information 
on molecular relaxation through a study of microwave 
absorption. Experimental results will be compared with 
the results of computations to evaluate proposed models 
of relaxation. 


experimental data, the model can be proved or dis- 
proved—an aid in establishing theory—and relaxation 
times can be determined. 

The behavior of molecules and their interactions is 
of considerable interest in basic research and in vari- 
ous electrical and electronic applications. One way 
of obtaining information on molecules is to study mi- 
crowave absorption in compressed gases containing the 
particular type of molecules in question. The mathe- 
matical formulation of one of the models under study 
was proposed by A. A. Maryott and E. Lustig, of the 
chemistry laboratory, in terms of a set of differential 
equations. The determinantial equations resulting 
from the corresponding matrices can be expressed in 
trigonometric form, and their solutions, necessary for 
the explicit line shape expression, can thus be readily 
obtained. 


Molecular Collision Studies 


The Heat Division is interested in the probability of 
transfer of translational to vibrational energy in col- 
lisions between atoms and diatomic molecules. Re- 
cently a model was proposed by K. E. Shuler and R. W. 
Zwanzig which assumed a hard-sphere interaction po- 
tential between the atom and the diatomic molecule in 
a collision where the atom approaches along the axis 
of the molecule. This problem was recently pro- 
eramed and carried out by R. Zucker of the Computa- 
tion Laboratory. 

The problem is formulated as an infinite set of 
linear equations with infinitely many unknowns, whose 
coefficients are given as polynomials of increasing 
order. Of course, even an electronic computer has to 
approximate such an infinite system by a finite number 
of equations—in this case, about 10 or 15 were neces- 
sary. The unknowns solved for are transition proba- 
bilities from a given vibrational state to other vibra- 
tional states, and the computing machine is made to 
stop when the sum of the computed probabilities 
converges to unity. 


The solution depends on the initial state of the 
diatomic molecule, on its initial vibrational energy 
relative to the translation energy of the single atom, 
and on the ratio of their masses. Solutions for several 
combinations of these parameters have been obtained. 


Reactor Design 


To keep reactor design abreast of technological ad- 
vances, a careful analysis of the reactor core, where 
fissionable material is placed, is necessary. Entering 
into this analysis are equations which describe neutron 
behavior, and whose solutions give the neutron flux dis- 
tributions throughout the reactor and the amount of 
fissionable material needed to sustain the nuclear chain 
reaction. It is possible to solve these important equa- 
tions in one dimension by hand if the geometry is 
simple, but this is very time consuming. Furthermore, 
the one-dimensional solutions do not represent the 
situation well enough to establish the necessary design 
criteria. For these reasons, automatic computing 
machines are used extensively in designing nuclear 
reactors. 

The Bureau, now in the preliminary stages of design- 
ing a high-flux research reactor, is using high-speed 
computation in the core design. Of the many codes 


available for solving core equations, two are being em- 
ployed. One solves the equations in one dimension and 
is called “WANDA,” + and the other is a two-dimen- 
sional code called “PDQ.” ? It would take 2 or more 
days to solve the one-dimensional problem by hand 
whereas WANDA requires only 2 minutes on the high- 
speed computer. Such speed not only provides overall 
economy but also permits a much greater variation of 
conditions to be studied. PDQ requires about half an 
hour of machine time, but solves the two-dimensional 
problem which would be impossible to solve by hand. 
PDQ is being used to determine the critical mass of 
fissionable material and the neutron fluxes given by 
different geometrical configurations. WANDA is used 
to study the effect of small variations from the basic 
geometries established with PDQ. This work is being 
carried out by R. S. Carter. 


* WANDA, a one-dimensional few group diffusion equa- 
tion code for the IBM-704, by O. J. Marlowe, C. P. Saal- 
bach, L. M. Culpepper, and D. S. McCarty, Westinghouse 
Atomic Power Division TM-28 (November 1956). 

* PDQ, an IBM-704 code to solve the two-dimensional 
few-group neutron-diffusion equations, by G. G. Bilodeau, 
W. R. Cadwell, J. P. Dorsey, J. C. Faircy, and R. S. 
Varga, Westinghouse Atomic Power Division TM-70 
(August 1957). 


Automatic Programing Studies 


4 ie INSURE the effective operation of its own com- 
puting facility and provide basic data of use to the 
entire field of automatic computation, the Bureau is 
engaged in a series of studies on automatic program- 
ing. This work, under the direction of J. H. Wegstein 
of the Computation Laboratory, is aimed at developing 
techniques for the more efficient communication of 
problems from the formulator—a research scientist or 
an administrator for example—to a computer. 

An essential step in this direction is reducing the 
time required to program and code-check problems. 
The Bureau therefore concentrates on this point using 
three main approaches: It shares programs with other 
computer users, eliminating the time and effort re- 
quired in writing a new program; it applies user- 
oriented methods, which, in effect, incorporate much 
of the detailed coding into the original statement of the 
problem; and it works on automatic operation systems, 
which reduce manual computer operations by means 
of a supervisory code. 

Advancements along these lines are immediately put 
to use in the Computation Laboratory, where scientific 
and data processing problems are solved on a high- 
speed computer both for the Bureau staff and for other 
Government agencies. On a broader scale, results are 
made available to other laboratories through national 
and international cooperative activities and consulting 
services. 


Program Sharing 
Work by Bureau programers in automatic program: 


ing has been greatly reduced by a cooperative organ- 
ization called SHARE in which computer users ex- 
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change subroutines and other computer programs. 
The Bureau has contributed several programs to the 
SHARE library and has in return made use of pro- 
grams written in other laboratories. 

The success in sharing programs in computer-user 
organizations has pointed up the desirability of uni- 
versal sharing which would be independent of com- 
puter type. However, in order to do this, a common 
language is required. The Bureau is therefore en- 
gaged in international efforts to define common ma- 
chine languages in both the scientific and business 
data processing fields. 


User-Oriented Efforts 


A relatively new area of automatic programing de- 
scribed as user-oriented efforts is being effectively 
applied at the Bureau. Instead of writing problems in 
the difficult language of the computer, programers now 
write in a pseudo language,' which the computer trans- 
lates into its own language. ‘This allows large prob- 
lems to be coded more easily and makes practical the 
automatic solving of small “one-shot” problems. 
However, one of the obstacles to the efficient use of a 
large computer on small problems is the communica- 
tion between problem author and the programer. AI- 
though the author may have a problem well organized, 
he may not express it in a form that is convenient for 
computer programing. 

To remedy this situation the Bureau has established 
a service called the “Black Box Computer” designed 
by J. H. Wegstein and G. M. Galler of the Computation 
Laboratory. In this operation, the computer, its 
operators, card handlers, key punchers, and coders 
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The “Black Box Computer”—a tool devised to improve communication between problem formulators and the auto- 


matic computer. 


Ruth Davenport of the dental research laboratory submits her problem, written in such a way that 


a programer’s services are not required, to S. Prusch, who writes the necessary portions in FORTRAN machine 


language. 


them to the computer with G. M. Galler at the console. 


operated by P. O’Hara, and returned to the scientist. 


The problem then goes to card-punch operators Betty Pailen (left) and Mary Schultz (right) and from 
Finally, results are printed out on an offline printer, 


make up the “computer”. The service is effectively a 
tool for the formulator of the problem who states his 
problem almost as he would naturally but in such a 
way that the services of an expert programer are not 
required to translate it into conventional computer lan- 
guage. Instead a coder takes the problem in its 
original form and utilizes previously programed pack- 
ages as well as other automatic programing aids in 
handling it. 

The most useful program in the Black Box Computer 
so far is called Tablemaker. It is a general purpose 
program for reading, computing, and printing tables 
of numbers. In addition to computing tables of arbi- 
trary functions, this program can numerically inte- 
grate, differentiate, interpolate, and difference. One 
of its outstanding features is its ability to fit curves 
when a table of numbers X and a table of numbers Y 
are given, representing points on the X—Y plane. The 
author also provides the function to be fitted in the 
form Y=aj,(x) +bf2(~) +cfs(x) + ..., where all 
the functions of X are arbitrary functions, and a, b, 
c, .. . are the parameters to be determined. Another 
useful feature of the program is its curve plotting 
routine. Tables of X and Y representing a curve in 
the X—Y plane are read by the computer and the curve 
is then plotted by characters printed on the off-line 
printer. 


Automatic Operator Systems 


Another area of automatic programing receiving 
emphasis is the use of automatic operator systems. In 
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addition to input translaters, compilers, libraries of 
subroutines, post-mortem routines, and output con- 
version routines, such a system provides a supervisory 
code to minimize the amount of “button pushing” re- 
quired of computer operators. This code automati- 
cally provides programers with information that would 
otherwise have to be obtained manually. For example, 
if the computer stalls in a program, the supervisory 
code can often provide enough information for the 
programer to determine the cause of the trouble. 

An automatic operator system, called Corbie, was 
developed by J. H. Wegstein and J. Cooper in 1957. 
This system is still in use and has provided much 
valuable help to programers; it has also been a valuable 
research and experimental tool in automatic program- 
ing. Unfortunately the Corbie system was not suited 
to problems run on magnetic tape. During the past 
year (1959) the computer was expanded to include 
32,000 words of high speed memory and an off-line 
machine for printing from magnetic tape, making it 
necessary to use a much more elaborate automatic 
operator system. With only minor changes, the Bell 
Laboratories’ system was adapted by G. M. Galler for 
use at the Bureau. With further additions, it is ex- 
pected to supersede the Corbie system entirely. 


*An example of this is the automatic program called 
FORTRAN, for Formula Translation. For further in- 
formation on computer language, see From formulas to 
computer oriented language, by J. H. Wegstein, Commun. 


of the ACM, 2, 6 (1959). 


OMNIFORM I—A General-Purpose Program for High-Speed Computers 


spe Bureau has developed a general-purpose pro- 
gram which enables scientists to apply high-speed 
computation techniques directly to many of the prob- 
lems arising in research. This program called OMNI- 
FORM I, was designed jointly by G. M. Galler of the 
Computation Laboratory and J. Hilsenrath of the 
thermodynamics group to eliminate some of the delay 
caused by a wide variety of tedious computations which, 
because of their ad hoc character, are not usually 
programed. 

OMNIFORM I is an interpretive program which en- 
ables even a user unfamiliar with computer programing 
techniques to instruct the machine to perform special- 
ized calculations. Although this program was created 
primarily to meet the needs of workers in thermody- 
namics and molecular physics, the mathematical 
operations involved are applicable to calculations in 
many fields of physics, chemistry, and engineering 
where properties of materials or other tabular data are 
required as a function of one or two variables. 

The program has anticipated a wide variety of 
algebraic manipulations of the type performed in many 
laboratories at the Bureau and elsewhere. The func- 
tion types built into the program are identified and 
called for by numbers similar to the stock numbers used 
in shoe stores to identify shoe styles. Thus, if a scientist 
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Calculations from one Bureau laboratory which were per- 
formed on OMNIFORM. Individual programs written 


forsee calculations would total approximately 2,500 
cards. 
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needs to substitute a set of members in a power series, 
he calls for a function of type 1—a power series of the 


form F(1)=>‘a,X°s: tells the machine how many 
j=1 


terms, n; and supplies n number pairs, a;, b;. He 
also supplies the arguments desired, information con- 
cerning further manipulations, if any, on the functions 
generated, and instructions for the print format. 

One manipulation which can be performed permits 
the taking of powers, roots, logarithms, antilogs, and 
exponentials of the computed function values. These 
operations are performed by an equation called a, 
“Transformer”, for which instructions are supplied by 
a set of seven numbers, usually integers. The results 
so modified may be operated on further in a “Compiler” 
and an “Integrator”. 

Additional features of the program permit a number 
of transformations of the variable; read-in of tabular 
values for combination with computed results; cascad- 
ing of operations such as is involved in continued frac- 
tions or roots-within-roots; choice of print format (any 
number of columns up to seven may be had in either 
fixed- or floating-decimal form) ; and optional print- 
out of results at various stages during computation. 

At present the program has built into it 16 function 
types: power series—with positive, negative, or frac- 
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tional powers; natural and common logarithms; ex- 
n 

ponential terms such as Sha; Xiexp (eX) ; 6 
ik 


ja 
trigonometric functions; and 6 hyperbolic functions 
(table). Other special functions which are contem- 
plated for future inclusion are: Bessel functions, and 
Legendre. Laguerre, Hermite, and Tschebyscheff 
polynomials. 

The program was assembled using the FORTRAN 
programing system with subprograms from the FOR- 
TRAN Library of routines. In addition, specialized 
subprograms were developed to provide a complete 
program giving utmost flexibility in input, computing, 
and output capabilities. These programs are a perma- 
nent part of OMNIFORM and hence provide a self- 


contained computing deck suitable for immediate use. 


The deck, comprising about 500 binary punched cards, 
can be read by the computer in 2 min, after which the 
computer is ready to interpret the 10 to 30 cards of 
input data which direct it to the individual problem. 

OMNIFORM is primarily a table generator. It can 
produce a table of up to 50 entries (50 different func- 
tions) for each of 100 arguments. Among the tables 
which have been computed on OMNIFORM thus far 
are: Ideal gas thermal functions for neon and electrons, 
anharmonic corrections for diatomic molecules, vapor 
pressure of solid nitrogen, virial coefficients for nitrogen 
gas, negative exponentials corresponding to certain 
atomic energy levels, polynomial substitution, table 
conversion to other units, table editing (column 
switches, etc.), and various ad hoc mathematical tables 
for selected arguments. 


Handbook of Mathematical Functions 


A! the request of the National Science Foundation, 

the Bureau is assembling an extensive collection of 
mathematical tables.t_ This work, called the Hand- 
book of Mathematical Functions, is designed to cover 
almost all the known transcendental functions. It 
was begun by the late M. Abramowitz, Chief of the 
Computation Laboratory, and is now under the direc- 
tion of I. A. Stegun also of that group. In addition 
to numerical tables, the handbook will include graphs, 
formulas, and descriptive text. 

The need for such a handbook has by no means been 
eliminated by the development of the high-speed com- 
puter but has perhaps even been intensified. Tables 
are required not only by scientists and engineers who 
do not have access to machines, but also for prelimi- 
nary surveys of problems prior to programing for 
machine solution. 

The handbook designed to meet these needs is to 
be a 1200-page volume in which 700 pages will be text 
and graphs, and 500 pages will be tables. The book 
will cover trigonometric and hyperbolic functions; 
sine, cosine, and exponential integrals; gamma func- 
tions; the error function; Fresnel integrals; Legendre 
functions; Bessel functions and integrals; Struve func- 
tions; hypergeometric functions; elliptic functions and 
integrals; parabolic cylinder functions; Mathieu fune- 
tions; spheroidal wave functions; orthogonal poly- 
nomials; Bernoulli numbers; binomial coefficients; 
the Riemann zeta functions; combinatorial analysis; 
interpolation; differentiation and integration; prob- 
ability functions; miscellaneous functions; scales of 
notation; and Laplace transforms. Reference tables 
of the mathematical and physical constants and of 
powers and roots will also be included. 

The Handbook of Mathematical Functions will bring 
together for the first time material from many separate 
sources. For example, it will include tables published 
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in 54 volumes of the Bureau’s Applied Mathematics 
Series, 13 volumes of the Columbia University Press 
Series, and 37 volumes of the Mathematical Tables 
Series, as well as newly computed tables and sections 
of tables from other sources. Most of this tabular 
material will be given to 5 decimal places or 5 signifi- 
cant figures and, for the more basic functions, values 
from 8 to 25 places will be given. In most cases, the 
values of the function will cover the entire range of the 
argument. 

Material is being contributed by many experts in the 
field of special functions. The Bureau is represented 
by contributions from M. Abramowitz; P. Davis, Chief 
of the Bureau’s Numerical Analysis Laboratory; W. 
Gautschi, now with the Oak Ridge National Labora- 
tory; M. Newman and K. Goldberg, of the Numerical 
Analysis Laboratory; E. Haynesworth, of the Com- 
putation Laboratory; M. Zelen, of the Statistical Engi- 
neering Laboratory; and N. Severo, now at the Uni- 
versity of Buffalo. Contributors from other labora- 
tories include G. Blanch, of the Wright Air Develop- 
ment Center; U. Hochstrasser, of the Swiss Embassy; 
A. N. Lowan, of Yeshiva University; Y. Luke, of the 
Midwest Research Institute; J. C. P. Miller and L. 
Slater, of Cambridge University; L. M. Milne-Thom- 
son, of the Mathematical Research Center; F. Oberhet- 
tinger, of Oregon State College; F. Oliver, of the Na- 
tional Physical Laboratory; and T. Southard, of the 
University of California at Los Angeles. _ Textual ma- 
terial will treat important properties of functions in- 
cluding definite and indefinite integrals, infinite series, 
asymptotic expansions, inequalities, continued frac- 
tions, differential equations, recurrence relations, poly- 
nomial and rational approximations and interrelations. 


*Problem solving on the high-speed computer, NBS 
Tech. News Bul. 42, 71 (1958). 


HIGH-SPEED COMPUTERS 
IN 


RADIO PROPAGATION 
RESEARCH 


A STORED-PROGRAM commercial digital com- 
puter is helping the Bureau carry out its radio 
propagation research program at the Boulder (Colo.) 
Laboratories with increasing speed, efficiency, and 
thoroughness. As knowledge of computational tech- 
niques advances, the computer is becoming a more 
powerful tool for data reduction, data analysis, and 
data-to-theory comparison. Problems in radio propa- 
gation illustrate one of the large number of fields in 
which the potentialities of digital machines are being 
explored. 

The growing indispensability of the high-speed com- 
mercial machine is clearly shown by the variety, extent, 
and importance of the problems undertaken at Boulder. 
Complex calculations, extensions of theories, and the 
production of essential graphs and tables are made 
possible or greatly accelerated by the use of a computer 
system. Equally important are the new lines of in- 
quiry or revisions of current theory that are frequently 
suggested by the results of a completed program. 

Solutions to many problems where the number of 
arithmetical operations may run into the millions are 
now almost routine. However, many others that re- 
quire more sophisticated computational techniques are 
being investigated. Many fundamental studies are 
being carried out with the aid of this computer, such as: 
Analysis of specular reflection from ionized meteor 
trails; worldwide prediction of transmission loss over 
ionospheric and tropospheric propagation paths; and 
the preparation of maps of ionospheric and tropo- 
spheric characteristics. 


Data Collection and Preparation 


Boulder Laboratories scientists are trying to develop 
recording techniques that will eliminate intermediate 
handling of data before processing by the computer. 
Unfortunately, present practical limitations preclude 
recording in form suitable for direct input to the com- 
puter. Intermediate steps of data preparation are re- 
quired; reduction or elimination of these intermediate 
steps is the subject of continuing effort. 

Reducing filmed records produced by ionospheric 
vertical sounding equipment is presently done by read- 
ing and scaling by the operator. The values are hand- 
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out of surface-wave diffraction cal- 
culations on a machine tabulator. A 
versatile adjunct to the commercial 
computer system, the tabulator or 
line-printer, depending on how it is 
wired, will put data into the system, 
print out results, and perform inde- 
pendent addition operations. 


punched on cards. These cards, in turn, are used as 
input data for the computer. An integrated procedure 
using standard format, a system of machine cards, and 
routine manipulation of data has improved the speed 
of this type of data reduction. 

At present, data taken in the field are recorded in 
several ways beside the familiar strip-chart techniques. 
These include punched cards or punched paper tape 
for digitalized data, and 12-channel magnetic tape 
analog recordings. The analog data are either fed 
directly into special analog computers or later con- 
verted to the digital form required by the computer. 
Recording data from strip-charts to punched cards can 
be accelerated by using a machine on which the 
operator alines crosshairs, pushes a button, and the 
machine records the data in card form suited for com- 
puter use. 

Most of the data in radio propagation experiments 
are continuously recorded over a period of time; rela- 
tively few data are taken on a short-term sampling 
basis. Moreover, most of the information is in analog 
form instead of the digital form required by the present 
computer. Thus, where a large amount of data reduc- 
tion is required, to not more than between 1 and 5 
percent accuracy, special purpose analog computers 
have been constructed. One example is the special 
application of an automatic wave analyzer at Boulder 
Laboratories. Compatible with the 12-channel magnetic 
tape used to record the data, this computer plays back 
the tape and produces a spectrum analysis of the wave 
function. Other analog computers have been developed 
or modified to analyze amplitude distribution, fade rate, 
and fade-rate duration. High-speed digital voltmeters 
are used to convert this analog information to digital 
form for use in the central computer. 


Data Analysis 


Analysis of data by computer systems is of particular 
value to some radio propagation research studies since 
it permits meaningful summaries of large amounts of 
statistical data. An example is the development of the 
NBS standard radio refractivity index, recommended 
for international adoption by the International Consult- 
ative Committee on Radio (CCIR). A quantitative 
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W. V. Mansfield watches the print- | 


Electroplotter, being set up by Maurine Smith, expands 
usual computer system operations to obtain graphic rep- 
resentations of various types of data. Information is fed 
into the electroplotter by cable connected to any one of 
three machine-room components of the computer system: 
Card punch, summary card punch, or tabulator. 


understanding of the bending of radio waves is impor- 
tant to the successful use of radar height-finding 
methods, tropospheric forward-scatter predictions, and 
radio-guidance systems. 

Similarly, analysis and comparison of data by com- 
puter are assisting materially in studies of an experi- 
mental radio communication system that will use 
reflections from ionized meteor trails for two-way 
message transmission. Early results indicate that inter- 
mittent meteor burst communication can compete effec- 
tively with other long-distance systems, and that it is 
relatively free from ionospheric disturbances. 

Another data-to-theory comparison in the field of 
sferics (naturally occurring atmospheric radio signals) 
involved the analysis of these signals into individual 
frequency components. This analysis proved valuable 
in studies of the attenuation of sferics and low-fre- 
quency radio waves traveling through the atmosphere. 


Improvement in Prediction Services 


One of the Bureau’s services is the prediction, several 
months in advance, of ionospheric characteristics which 
permit the selection of appropriate transmission fre- 
quencies for long-distance propagation paths. The 
computer’s ability to perform large amounts of cal- 
culation involved in the theoretical models allows a 
testing of these models against observed data on a scale 
which would not be possible using hand or desk 
calculator techniques. 

Present prediction services for ionospheric radio 
propagation techniques are undergoing a change. In 
the future, MUF (maximum usable frequency) predic- 
tions are expected to issue in two forms: Printed graphs, 
as at present, for users without computers; and on mag- 
netic tape or on punched cards for those who do have 
computing facilities. It is proposed that users of the 
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latter form of prediction service feed this packaged 
data into their computer along with the coordinates of 
the paths to be used. Their computer would then 
develop a detailed prediction for these paths. At 
present, the prediction process requires a hand operation 
in reading values from the prediction graphs. 

The Bureau has developed a computer program 
which can be used to produce contour maps of the 
ionosphere. The larger computer at the Washington 
Laboratories had to be used even in the early stages 
since the one available at Boulder did not have suft- 
cient capacity for the task. These maps are based on 
median values of foF2 for 24-hr periods, tabulated for 
each station in its own time zone during a given month. 
It is expected that eventually similar maps can be pre- 
pared from other ionospheric characteristics, meteor- 
ological data, or other geophysical quantities. 

In this task of producing numerical ionospheric 
maps, data are analyzed by progressively fitting series 
of orthogonal functions in three coordinates: time, 
latitude, and longitude. The process requires smooth- 
ing for the clustering of points and noise interference. 
When completed, the time and geographical variations 
of foF2 are represented by a table of coefficients for a 
quite complicated function of three variables, equiva- 
lent to a “mathematical map”. 

This project, however, is in no sense a mechanization 
of the Bureau’s present prediction program. It is 
instead an attempt to make the most effective use of 
high-speed computers in answering certain propagation 
problems. As a byproduct, this type of computer pro- 
gram offers a good tool for further ionospheric 
research. 

Substantial changes in the scope and accuracy of 
predictions for other types and ranges of radio propa- 
gation are also in process. Those concerned with the 
development of VHF and UHF equipment (30 to 3,000 
Mc), will soon have the benefit of tables of hourly 
median levels of transmission loss over predetermined 
paths. The computer has made possible the considera- 
tion of models involving many more parameters. Pro- 
vision within the model can now be made for air 


Louise Kiniston keeps track of the operation of auxiliary 
components of the computer system: A single-column 
automatic punched-card sorter (left); and a statistical 
sorter, which can sort many columns at the same time 


(right). 
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turbulence and other weather conditions. The data 
are based on a 5-yr observation program over 300 main 
paths in the United States. Variations of the data 
from previous theoretical models are compelling re- 
consideration and correction of these models. 


Tables and Graphs 


Production of tables and graphs of special functions 
is another phase of computer work at the Boulder 
Laboratories. The ability of the computer to handle 
very complex mathematical operations allows the de- 
velopment of entirely new functions, as well as modi- 
fications of those already known for particular research 
or general use. 

Finally, many theoretical calculations have been ex- 
tended far beyond the limits of previous work by using 
the computer system. Extended theoretical calcula- 
tions and the machine plotting of graphs were used in 
a pilot study of the practical value of an atlas of dif- 
fraction curves, based on the Van Der Pol-Bremmer 
theory of surface-wave diffraction. A larger com- 
puter is indicated to complete a large-scale compilation. 
These curves depend on a variety of parameters and 
can be used to give the field strength that can be ex- 
pected from transmissions at frequencies from 10 ke 
to 10,000 Mc, using the diffraction method of radio 
propagation. 


During an analysis of ionospheric radio propagation data 
by digital computer, Bonnie Laubach (left) checks the 


operation of the punched paper tape reader. This unit, 
designed and built to the Bureau’s specifications, con- 
verts information carried by tape—up to 7 bits wide— 
at the rate of 150 characters per second into form the 
computer can use. Patricia Lollar (right) monitors the 
control console and the unit which houses the magnetic 
storage drum, 
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Amplitude Distribution Analyzer System, one of several 
special-purpose data analyzers adapted or modified from 
commercial equipment at the Boulder Laboratories. 


Scheduling of Equipment 


Currently, computer equipment at the Boulder Lab- 
oratories is in use for three shifts, 5 days a week, with 
some weekend use. Eight hours weekly for preventive 
maintenance, inspection, and troubleshooting provides 
a break in scheduled programing. Work of an identi- 
fied trouble area, on an unscheduled basis, can often 
be performed while data are being processed in other 
units of the system. A test set of cards provides a 
complete check of every unit in the computer group 
in half an hour; a test paper tape checks data handling 
from 7-bit binary code tape. 

Two methods of operating the computing equipment 
are in regular use. The computing room staff regu- 
larly rus the longer and more routine problems dur- 
ing the evening and night hours. The daytime hours 
are available to individual scientists to check out new 
programs or make special runs where they may feel 
that they can contribute personally. Individual famil- 
iarity with computer operations appears to have a 
special value in that it qualifies the scientist to plan 
integrated programs for reduction, analysis, and pro- 
jection of data, rather than working out the program 
piecemeal. 
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Frequency Selection for Aurlane 


“Road Markers’ a 


4 es Bureau, at the request of the Federal Aviation 
Agency, is studying automatic computation methods 
for determining the best possible operating frequen- 
cies' for radio transmitters used as “road markers” 
on airlanes. 

A network of such transmitters marks routes between 
cities in the United States by sending out signals which 
assist pilots in flying a straight-line course. The 
rapid expansion of commercial and military air opera- 
tions makes it necessary to add a substantial number of 
new transmitters to this network each year. The 
locations of these transmitters are determined by tech- 
nical and economic considerations. However, the 
choice of carrier frequencies constitutes a surprisingly 
difficult combinatorial problem. Efforts to solve this 
problem are being carried out by L. Joel, G. M. Galler, 
and A. J. Goldman of the Applied Mathematics 
Division. 

The difficulties of frequency selection stem from the 
fact that transmitters with identical or neighboring 
carrier frequencies must be spaced widely enough to 
prevent signal interference. Furthermore, this must 
be accomplished within the range of 100 discrete fre- 
quencies assigned to FAA. However, not only is the 
“interference radius” of a transmitter large (approxi- 
mately three times the radius its signal “serves”) but 
the number of transmitters in existence is already con- 
siderable and is increasing rapidly. For these reasons, 
assigning a frequency to a new transmitter without 
introducing any interference requires a laborious ex- 
amination of many “old” transmitters. 

Such an assignment may in fact be impossible with- 
out changing the frequency of one or more old trans- 
mitters. This, in turn, may create new interferences 
in the system and require alteration of the frequencies 
of still more transmitters. The insertion of a new 
transmitter into the network has sometimes required 
frequency changes at as many as 11 old transmitters. 
Such changes are expensive and disturb the smooth 
operation of the system, since pilots must be informed 
of and become accustomed to the alterations. An addi- 
tional requirement introduced into frequency assign- 
ment is therefore the limitation of the number of 
changes. 

As the frequencies and transmitters (old and new) 
are finite in number, there are only a finite number of 
ways in which frequencies can be assigned throughout 
the network. This means that in principle the problem 
could be solved by examining all such network-wide 
assignments, rejecting those which lead to interference, 
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Typical station for broadcasting direction and distance 


information to pilots. The Bureau is developing auto- 
matic frequency selection techniques for the many new 
installations which the Federal Aviation Agency erects 
each year. Photo courtesy of FAA. 


and then selecting from the remainder, one which as- 
signs the original frequencies to the largest number 
of old transmitters. This procedure is impractical be- 
cause the number of cases to be examined, though 
finite, is so enormous (exceeding 10*°°°) that the in- 
vestigation could not be carried out in a reasonable 
time even on the most rapid electronic computer. 
What is needed is a systematic (rather than an ex- 
haustive) procedure. 

The Bureau’s efforts toward such a procedure can 
be roughly divided into two categories, which may be 
called ad hoc methods and model construction. The 
preliminary direction of the ad hoc activities was sug- 
gested by the fact that many new transmitters require 
changing at most one old transmitter. As a temporary 
aid to the FAA and an introduction to the general prob- 
lem, a routine has been written to make such “easy” 
frequency assignments when possible. 

The major activity in the model construction area 
has been the examination of a formulation which re- 
duces the frequency assignment problem to the maximi- 
zation of asum. Computational methods and computer 
codes, although not immediately applicable in this 
instance, are available for this general type of problem, 
and efforts are being made to adapt them for frequency 
selection. 


* “Best” in the sense described below. 
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MONG the first areas in which the potentialities of 
automatic computation were exploited and its 
techniques adopted on a broad basis was the solving of 
the large accounting, statistical, operational, and other 
specialized problems that arise in business, industry, 
and government. ‘The intense activity of this kind 
which still exists is immediately evident in the number 
of data processing and operations analysis problems 
the Bureau is called upon to solve. The demand for 
assistance in these relatively new areas is so great, that 
at present such work accounts for about 20 percent of 
the time the Bureau’s 704 computer is in use and about 
one-third of the coding effort expended. 

Data processing and operations analysis problems 
come not only from the Bureau’s staff but from other 
government agencies who call upon the NBS computing 
facility for help in applying present-day electronic 
methods to information processing, information re- 
trieval, accounting, reporting, and large-scale analysis. 
To illustrate the scope of this work, some of the prob- 
lems which have recently been solved or on which the 
Bureau is now working are discussed below. 


Autotag 


A war-game study called Autotag, which simulates 
battles involving tanks and antitank systems, was pro- 
gramed and carried out by the Bureau for the Combat 
Operations Research Group of Technical Operations, 
Inc. Use of computers to play war games in this way 
is becoming more and more frequent because the pres- 
ent development of weapons systems is so rapid that 
it is difficult to obtain the experience needed for compe- 
tent tactical decisions. The computer, with its ability 
to simulate hundreds of maneuvers cheaply and rapidly, 
can in many cases provide an acceptable substitute for 
actual military exercises. In addition, the computer, 
instead of acting intuitively, can analyze numerous 
possible courses of action and supply the proper one 
for the battle situation at hand. 

In the Autotag study, conducted by W. G. Hall and 
J. M. Beiman of the Computation Laboratory, a combat- 
simulation model was designed to evaluate the perform- 
ance of antiarmor weapons systems in a variety of 
topographic and environmental situations. In any one 
simulation, each side may have up to 31 weapons which 
may be any or all of 10 different types. Each weapon 
is specified by 69 adjustable parameters. Some of these 
specify weapon type. They include, for example, maxi- 
mum and minimum range, hit probabilities, probabi- 
listic delay times for weapon action, dispersion of fire 
as a function of range, and kill-if-hit probabilities. 
Other parameters relate to a particular weapon. Posi- 
tion and velocity coordinates, status indications (“dead” 
or “‘visible to enemy”, for example), and indices 
specifying the weapon’s target selection rules are some 
of these. 

Each side is also permitted up to seven mine fields 
and artillery fire at any of seven intensities before and 
during the battle. Two additional parameters specify 
the vulnerability of each weapon to mines and artillery. 
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It is assumed that no weapon is ever “killed” by the 
artillery or mines of its own side. 

Simulations are carried out in sequences which in- 
clude from 1 to 15 games, differing only in the random 
numbers controlling the outcomes of probabilistic 
events. The results are recorded on magnetic tape and 
an additional computer routine later performs the neces- 
sary tabulations and statistical analyses. Detailed 
histories of each simulation are also preserved for 
further analysis. 


Electrocardiograph Analysis 


Under the sponsorship of the Veterans’ Administra- 
tion, the Bureau is designing computer programs that 
will analyze human heartbeats. The heartbeat data 
are in the form of electrocardiograms which have been 
converted to digital readings by a method developed 
by the Bureau’s Data Processing Systems Laboratory." 
The ability to provide a rapid, detailed analysis of 
such records suggests a method for improving the pro- 
cedures employed in diagnosing various heart ailments. 

The programing of the computer for this purpose is 
being carried out by R. J. Arms of the Computation 
Laboratory. His immediate objectives are to detect 
the time intervals of important phases of the heartbeat 
cycle, to locate the least-squares planes in these cycles 
and determine the goodness of plane fit, and to find 
mean vectors, mean velocities, and the velocity vari- 
ance. Although the time-interval detection has been 
found difficult and may not be suitable to digital 
methods, finding the means and the least-squares fits is 
a common digital problem. Digital methods are also 
immediately applicable to the analysis of these param- 
eters along with about 100 others, which will include 
age, weight, and known diseases, to name a few. 


Low-Rent Housing Analysis 


Under the sponsorship of the Public Housing Ad- 
ministration, the Bureau analyzes for error the yearly 
reports on tenants in federally aided or owned low-rent 
housing projects. This annual status examination is 
necessary to determine the tenants’ eligibility for con- 
tinued occupancy. Approximately 400,000 such ex- 
aminations must be centrally edited for inconsistencies 
and prepared for statistical presentation. 

Inconsistencies and errors may occur among almost 
any combinations of the many items of data required 
for each tenant. These data include four items on in- 
come, nine on exemptions, two on rent, three on family 
composition, and others on age, race, and service 
status, plus pertinent remarks where amplification is 
necessary. Also, the diversity of backgrounds among 
the individual respondents makes for further difficulty 
in predicting the volume and nature of errors likely to 
occur. It is therefore necessary to make numerous 
checks on the relationships among the several items of 
data. 

To facilitate these editorial operations, G. H. Urban 
and A. T. Nelms of the Data Processing Systems Lab- 
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oratory developed a comprehensive computer program. 
The initial step in preparing reports for input to the 
computer is the coding of remarks—a clerical opera- 
tion performed in accordance with a simple two-digit 
code. Data cards are then punched from the reports 
and converted to magnetic tape. An additional input 
to the program is the “project header file” which con- 
tains in one record per housing project all the rules 
eoverning the administration of that project, the name 
and the full address of the respondent, and various 
other project descriptions for statistical purposes. 

The main editing program contains about 5,800 com- 
puter instructions. Here the data are analyzed criti- 
cally for reporting inconsistencies and policy viola- 
tions. Whenever an error is detected, it is recorded 
on magnetic tape. Six machine words are required to 
describe and index the error. Whenever possible, the 
program corrects erroneous data and derives entries 
for items left blank. The data are further distributed 
into meaningful statistical arrays. 

Another program sorts the errors within each project 
by the assigned query numbers and writes an exposi- 
tory paragraph describing the inconsistency or viola- 
tion. These paragraphs become the substance of 
letters automatically generated and addressed to ap- 
propriate regional and project authorities. The names 
of the families involved are tabulated for printing fol- 
lowing each paragraph. ‘The letters are transferred 
to magnetic tape and are then ready for off-line pro- 
ductive printing on continuous-feed letterhead paper. 


Payroll Operations 


A program that will allow the Bureau to use a com- 
puter for generating payroll and other personnel serv- 
ices reports has recently been completed by M. L. 
Paulsen and P. L. Ruttenberg of the Computation Lab- 
oratory. This development was based on information 
obtained in an original study which G. Urban of the 
Data Processing Systems Laboratory conducted under 
the Bureau’s program of assistance to other Govern- 
ment agencies. The complexity of the present under- 
taking is indicated by the many reports and documents 
which must be considered. These include the master 
pay list, the payroll for personnel services, check issues 
and statements of earning, timecards, savings bond 
schedules and purchase lists, earning reports, and with- 
holding tax statements. 

Like many data processing problems, the payroll 
operation utilizes a master tape, which in this case con- 
tains, for each employee, all pertinent information 
(name, pay rate, number of dependents, etc.). This 
master tape is continually updated by merging it with 
records of personnel actions such as new appointments, 
promotions, and separations. From the updated 
master tape and the timecard information for a given 
pay period, the computer produces payroll reports. 

Test runs have shown that the computer can dupli- 
cate in a few minutes the work which it now takes ap- 
proximately 80 man-hours to accomplish each pay 
period. 


16 


oe Controid 


A program which allows a computer to print out “trees” 
of minimum traveltime paths is used to make runs for 
the District of Columbia Department of Highways and 
Traffic. The tree shown illustrates how traffic behaves in 
approaching zone 36 in the heart of the Washington 
(D.C.) business district. 


Traffic Studies 


The Bureau is engaged in programing and solving 
several data processing problems on traffic for the 
Bureau of Public Roads.2 This work, carried out by 
S. T. Peavy and R. N. Varner of the Computation Lab- 
oratory, is basic to the effective design of new roads 
and highway networks. 

One important aspect of the investigations is the 
development of equations that will accurately predict 
trafic patterns. The Bureau is presently analyzing 
data on traffic passing through a signalized intersection 
in efforts to derive such equations. Part of the coding 
for this problem and an auxiliary code for the analysis 
of data on highway accidents have been completed. 
The results of the latter code, which relates the accident 
rate to control of road access, were requested by the 
Congress for a highway safety study it has authorized. 

Another problem, dealing with the assignment of 
projected traffic to existing and proposed routes, is 
handled by a program which computes a “tree” of 
minimum travel-time paths from each point in an urban 
area to all other points. Bureau facilities are used to 
make runs for the Washington (D.C.) area, the 
Minneapolis-St. Paul area, and the Columbus (Ohio) 
area. 


Linear Programing 


The Bureau solves a variety of logistics problems 
quickly and easily by making use of a relatively new 
mathematical field known as linear programing. In 
applying linear programing to a problem the first step 
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Ida Rhodes translates a Russian sentence into English. 
Mrs. Rhodes is developing automatic translation tech- 
niques which will allow a computer to make use of the 
syntactical relationships between words—thereby closely 
paralleling the thinking of a human translator. 


is to interpret mathematically the various conditions 
imposed. On the basis of this analysis, one of several 
linear programing computation procedures, previously 
coded for the computer, is applied to the data of the 
specific problem. 

The applicability of linear programing has greatly 
increased with the advent of high-speed computers, 
whose size and speed permit the treatment of larger 
and more complicated systems. Furthermore, the in- 
creasing availability of such advanced computers sug- 
gests that the utility of linear programing can be 
greatly extended. Already the Bureau is considering 
problems which are formally similar to those solved 
by linear programing techniques but have previously 
been thought to be beyond the scope of this method. 
Thus, Bureau mathematicians A. J. Goldman and L. S. 
Joel are investigating the reduction of very-large-scale 
systems by approximate and combinatory methods, 
nonlinear programing, and stochastic programing. 

The Bureau is attempting to extend linear program- 
ing to a variety of problems ranging from the selection 
of frequencies for air navigational radio stations at the 
request of the Federal Aviation Agency * to the alloca- 
tion of medical supplies to hospitals through a depot 
system with variable unit costs of shipment and stor- 
age—a problem submitted by the Veterans’ Admin- 
istration and the General Services Administration. 
Other work includes two studies for the Federal Com- 
munications Commission. In one, the Bureau is de- 
termining the optimum channel frequency and station 
location for commerical television stations. In the 
other, it is seeking the optimum trade-off time between 
service and interstation interference in AM and FM 
radio stations. 


Language Translation 


The Bureau is engaged in a program, sponsored by 
the Army Office of Ordnance, to develop mechanical 
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translation techniques.* The immediate objective of 
this work is to use an electronic digital computer to 
translate technical Russian into English. 

The mechanical translation studies, under the direc- 
tion of I. Rhodes of the applied mathematics staff, 
resulted in a process which goes much further than the 
word-to-word translation ordinarily associated with 
computers. This method is designed to parallel the 
thinking of a human translator. Furthermore, it ap- 
plies iteration to approach the correct syntactical inter- 
pretation of the Russian source sentence. The transla- 
tion proceeds in three main steps: First, relevant in- 
formation, stored in a machine glossary and auxiliary 
listings, is collected; second, source words, symbols, 
and punctuation marks are recognized by use of this 
stored information, and the clause and phrase structure 
of the source sentence is tentatively determined; third, 
a structural analysis of the source sentence with simul- 
taneous refinement of clause and phrase structure is 
made and the corresponding “pidgin” English is 
printed out. The heart of the method, which con- 
stitutes the unique contribution of the Bureau, is the 
prediction, after each source word, of future occur- 
rences on the basis of stored syntactical relationships 
between Russian words. 


Automatic Design of Computers 


The Bureau has devised techniques for using an 
existing computer to work out design details for new 
computers. The objective of this work, carried out 
under the direction of A. L. Leiner and A. Weinberger 
of the Data Processing Systems Laboratory was to 
mechanize some arduous manual design procedures. 

The procedure for designing a large-scale digital 
system involves four major steps: The preparation of 
system specifications, the organization of the system 
into large-scale functional blocks, the creation of the 
detailed logical design for each of these functional 
blocks, and the conversion of this logical design into 
detailed wiring plans. In deciding which steps of this 
procedure should be mechanized, the primary criterion 
was whether a machine or an experienced computer 
designer could turn out a better end product. In some 
cases it was found that a machine could not effectively 
compete with the empirical intuition of the experienced 
logical designer; in other cases, the use of a machine 
permitted a degree of design sophistication, accuracy, 
and efficiency not otherwise obtainable. On the basis 
of this criterion, the portion of the design procedure 
selected for mechanization was the conversion of a 
logical design into detailed wiring plans. 

The conversion process was carried out automatically 
for a new large-scale computer system. Suitable pro- 
erams were written by C. D. Coleman of the Data Proc- 
essing Systems Laboratory in collaboration with W. G. 
Hall, J. M. Beiman, and M. L. Paulson of the Computa- 
tion Laboratory. These programs provided a complete 
set of instructions for accomplishing the conversion of 
the computer logical design data. Methods for express- 
ing the logical design data in machine language were 


worked out by H. Loberman and W. W. Youden of the 
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The components of the heart vector, as seen on an oscil- 
loscope (left), are at present analyzed by 2-dimensional 
projections (right) of the generated space curve. The 
Bureau is in the early stages of devising digital tech- 
niques for carrying out the same analysis by machine. - 


blocks in the system, the computer automatically as- 
signed the numerous minor components to optimum 
physical locations with respect to the principal blocks. 
The computer then calculated the optimal wire routing 
between every connected set of components in the sys- 
tem, taking into account the restrictions on wire lengths 
imposed by their electronic and mechanical properties. 
Next, a topological diagram of the wiring tree, showing 
the physical location of each node and its connections 
with the other nodes in the tree, was printed out. 
Finally, the computer calculated and printed out the 
exact lengths of all wires, specifying terminal types, 
pin-to-pin connections, the text for labels on the wires, 
and a detailed inventory of the components and 
materials used. 

As an aid to maintaining the completed system, the 
machine also was used to generate detailed information 
on the interrelations between the logical structure of 
the system and the actual physical equipment used to 
implement the logic. The final output data consisted 
of 15,000 printed pages containing the complete wiring 
plans for the new system and source material for the 
maintenance manuals. 


Data Processing Systems Laboratory. The input data 
consisted of approximately 20,000 punched cards 
which the computer used to produce the complete wiring 
plans for the new system. 

The input transcriptions of the logical structure of 
the new system—expressed in terms of equations in 
Boolean algebra—were first checked rigorously for 
internal consistency. Then, using additional data de- 
scribing the physical locations of the principal building 


*High-speed computer to analyze human heartbeat, 
NBS Tech. News Bul. 43, 196 (1959). 

* Problem solving on the high-speed computer, VBS 
Tech. News Bul. 42, 71 (1958). 

’ Frequency selection for airlane “road markers”, VBS 
Tech. News Bul. See page 13 this issue. 

“Machine translation of Russian, VBS Tech. News Bul. 
43, 101 (1959). 


Automatic Computation in Problems of Number Theory 
coeficients of this series are precisely the desired 
numbers p(7). 
M. Newman of the numerical analysis group has 
made a study of the coefficients of the function 
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where the r’s are integers and the e’s are positive 
integers. A multilength code for the computation of 
these coefficients was written by E. Brauer and M. New- 
man and some 100,000 coefficients were computed. 
The results of these computations have been of ma- 
terial assistance in proving certain theorems in analytic 
number theory. For example it was proved that 
p[84n?—(n?-1) /24] is divisible by 13, provided the 
positive interger n is not divisible by 2 or by 3.1. In 
addition the computations verified and extended exist- 
ing tables of D. Ferguson and J. C. P. Miller (Cam- 


bridge, England) and suggested many new theorems 


LTHOUGH electronic computers are being used 
extensively to solve problems in science and busi- 
ness management, they are rarely applied to purely 
mathematical problems. The Bureau is one of the few 
organizations where such applications are being stud- 
ied, especially in the fields of number theory and com- 
binatorial analysis. A recent example of this work is a 
number theory study in which the coefficients of a 
“generating function” were computed. 

The method of generating functions is a classical 
procedure in number theory, in which a function is 
expanded into a power series—usually a routine 
process—and the resulting coefficients are used to solve 
a more difficult problem. For instance, let p(n) de- 
note the total number of partitions of n (so that 
p(5)=7 since 5 has the partitions 5, 4+1, 3+2, 
Saleh, ela Qari Dares eel iy ee eral) 


While it is difficult to compute p(n) directly for large 
n, one can derive a generating function from the “Euler 
product”, p(x) =(l-x) (l-«?) (1-28) ...; specifi- 
cally, if 1/p(x) is expanded into a power series, the 
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which were subsequently proved. 


* Further identities and congruences for the coefficients 
of modular form, by M. Newman, Canadian Journal of 
Mathematics 10, 577 (1958). 
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Target Selection by Aur 
Defense Systems 


A RESEARCH PROGRAM, sponsored by the Army 
Signal Corps, is being conducted by the Bureau 
to enable air defense weapons systems to realize their 
full potentialities. This work, directed by A. J. Gold- 
man, involves the development and implementation of 
rules for selecting the best possible target from among 
all the attackers in range. 

The more “glamorous” technological aspects of air 
defense systems sometimes tend to overshadow the fact 
that difficulties still remain after such a system is ready 
for use. For example, in the course of an air attack, 
the enemy may succeed in delaying recognition and 
accurate location of some attackers. Because such a 
“hostile environment” will prevent the defense weapons 
system from realizing its capabilities, it is extremely 
important to put each shot where it will be most effec- 
tive. The research of the Bureau’s applied mathe- 
matics laboratories is directed toward this goal. 

The factors involved in target selection include the 
nearness of an attacker to the defended area, and the 
velocity of an attacker. If all the relevant criteria 
point to the same attacker then of course there is no 
selection problem, but such a simple situation rarely 
arises. The function of a selection rule is to combine 
the relevant factors into a single composite criterion, 
which then yields an unambiguous choice of a best 
target among the attackers. Deciding precisely how 
the factors are to be combined (i.e., constructing a 
selection rule) requires a careful analysis of the rela- 
tionship between the various criteria (singly and col- 
lectively) andthe objectives of the defense. 

The task of target selection would be much simpler 
if each firing decision could be made with full fore- 
knowledge of the future course of the air defense battle. 
Exact knowledge of the future is obviously unattainable, 
while probabilistic prediction, though conceivably 
possible in principle, would require extensive calcula- 
tions quite incompatible with the requirement for rapid 
decision making. Because the future of the battle 
can be treated only in an approximate way, the effec- 
tiveness of a proposed selection rule cannot readily be 
determined by theoretical means. Experiments are 
therefore performed in which a high-speed computer 
mimics a number of air defense battles and reports the 
results of the target selections made. The Bureau was 


one of the first organizations to undertake such simu- 
lation studies. 


The Bureau is developing selection rules to enable air 
defense missile systems to choose the best target from 


among all attackers in range. The results of this work 
are applicable at installations such as the Army Air De- 
fense Command Post at Fort George G. Meade, Md. One 
part of this operation is the manning of a Missile Master 
Console (left), which controls the firing of missiles on 
friendly aircraft. Another console (center) supplies 
tactical direction, and a third (right) provides tactical 
monitoring. Photo courtesy of the Army. 


An examination has been made of the general prob- 
lem of conducting simulations as effectively as possible. 
The usual outputs of such studies are sufficient to indi- 
cate differences in effectiveness between two target 
selection rules, but cannot be used to detect the source 
of the difference. Methods for obtaining and present- 
ing this more detailed information (which is necessary 
for the systematic improvement of selection rules) have 
been developed and are currently being automated. 

Another study, which dealt with the computer func- 
tions common to most such simulation problems, sug- 
gested the creation of a computer program which would 
be “modular” in the sense that different selection rules, 
defense systems, and attack situations could readily be 
“plugged in” for investigation. A prototype* of such 
a program was prepared and is now in use at the 
Bureau. It performs the standard functions of simu- 
lator (“moving” the attackers and defense missiles 
along their paths), umpire (deciding “hit” or “miss” 
when a defense missile reaches its impact, and deter- 
mining when a target has carried out its mission), and 
historian (keeping a detailed record of all significant 
simulated events). The operation of this program has 
provided data on the possible use of computers to im- 
plement target selection rules. 


*FLEXSIM abstract, Flow charts and output format, 
distributed at the February 1959 meeting of the Com- 
puter Application to Military Problems group. Available 
from the CAMP representative at the Applied Mathe- 
matics Division of the National Bureau of Standards. 


Work on target selection rules will help air defense 
weapons systems, such as the four Nike-Ajax missiles 
shown, realize their full potentialities. Photo courtesy 
of the Army. 
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Statistical Planning of Physical Science Experiments 


for greater efficiency and precision 


4 Peees BUREAU conducts a continuing research pro- 
gram on methodology in the planning and analysis 
of physical science experiments. Studies have resulted 
in the development of several new families of experi- 
ment designs which provide maximum information 
with a minimum number of measurements. Work 
has also been carried out on the tabulation of experi- 
ment designs to allow even the nonspecialist to use these 
designs easily. 

The application of statistics and probability theory 
to the design of physical science experiments is a rela- 
tively new area in which the Bureau has had an active 
program during the past several years. The statistical 
design of experiments originally was applied to special 
problems in the agricultural and biological sciences. 
However, it became apparent that the use of sound 
statistical principles in the physical sciences could 
provide an improved basis from which to draw scien- 
tific conclusions. With the present emphasis on more 
precise and reliable methods of measurement, the 
interest in experiment design has been greatly 
intensified. 

Two types of designs receive attention at the Bureau: 
Fractional factorial designs, which are applicable to 
multifactor experiments; and block designs, which are 
employed where measurements form homogeneous sub- 
sets. The fractional factorial designs allow the scien- 
tist to plan his experiment so that it will include only 
a fraction of the possible experimental combinations, 
thus providing savings in time and effort. The block 
designs balance out the effects of systematic errors and 
relieve the experimenter of much of the burden of re- 
producing on every occasion exact environmental con- 
ditions and exact experimental units. 


Multifactor or Factorial Designs 


Experiments often involve several factors which 
jointly affect the material or process under investiga- 
tion. For example, in a chemistry experiment these 
factors may include temperature, time, concentration, 
and the presence or absence of a catalyst. Other 
things being equal, the smaller the number of factors, 
the fewer the difficulties encountered in the planning, 
analysis, and interpretation of the experiments. How- 
ever, because of possible complex interdependencies 
among the effects of the various factors, the problems 
of analysis and interpretation become difficult even for 
experiments with a moderate number of factors. 

Still another and perhaps more formidable problem 
is that as the number of factors increases, the number 
of different combinations of the factors increases ex- 
ponentially. Often the experimenter has neither the 
resources nor time to investigate all possible combina- 
tions of the factors. Instead he desires to select a 
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subset of all combinations which will serve as a basis 
for drawing conclusions about the effects of the various 
factors. 

The Bureau’s Statistical Engineering Laboratory has 
therefore tabulated a large number of fractional fac- 
torial experiment designs. Tables of fractional designs 
for experiments where each factor is varied at two 
conditions * or at three conditions * have been compiled 
with the partial support of the Army Chemical Corps 
and the Bureau of Ships, respectively. The latest work 
in this area has been the tabulation of fractional fac- 
torial designs for use when factors are varied both 
at two and three conditions. Unlike the other work, 
the theory for this mixed case was unknown. W. S. 
Connor, of the Statistical Engineering Laboratory, de- 
veloped an extensive mathematical theory for these 
“mixed” fractionals, which served as a basis for the 
actual construction of the designs.° 


Block Designs 


The statistical design of experiments was first applied 
in agriculture, where it was noticed that neighboring 
plots of land were more alike in their responses than 
plots that were widely separated. The principal of 
“blocking” was therefore developed in which experi- 
mental units which are close together are regarded as 
a “block.” Here “close together” may mean at the 
same time, on the same day, on the same materials, by 
the same operator, with the same apparatus, or what- 
ever is appropriate in a given situation. The principle 
behind the analysis is to provide the desired compari- 
sons using only the differences of measurements taken 
close together, that is, differences within a block. This 
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Fractional factorial experiment plan requiring only $l 
measurements out of a possible 243. Factors are indi- 
cated by letters; levels of factors by subscripts 0, 1, and 
2. The Bureau studies and develops such experiment 
designs to provide maximum information for a minimum 
number of measurements. 
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eliminates extraneous variations caused by differences 
between methods or processes. 

In agricultural field trials there is always some 
arbitrariness in the subdivision of a field into plots and 
the organization of these into blocks. However, in 
physical, chemical, and engineering experimentation, 
a natural block is usually indicated by the problem 
itself. For example, the measurements taken during 
a single day or the readings taken with a particular 
thermometer form a natural block. Furthermore, the 
number of experimental units in a block is often not 
arbitrary, but fixed by the nature of the problem. The 
two sides of a thermocouple and the six or eight cylin- 
ders of a car engine are illustrations. 

The Bureau’s continuing research on experiment 
design has provided several new classes of incomplete 
block designs which are especially suited for use in the 
physical sciences. Examples are the chain block design 
developed by W. J. Youden and W. S. Connor of the 
Statistical Engineering Laboratory; * the two-way elimi- 
nation chain block design of J. Mandel of the Organic 
and Fibrous Materials Laboratories;° and experiment 


designs with block size two, which are especially useful 
in high precision work.® 


*Fractional factorial experiment designs for factors at 
two levels VBS Appl. Math. Series 48 (1957). See also, 
Economy in the planning of experiments, NBS Tech. 
News Bul. 41, 78 (1957). 

* Fractional factorial experiment designs for factors at 
three levels, by W. S. Connor and Marvin Zelen, NBS 
Appl. Math. Series 54 (1959). 

* Construction of fractional factorial designs of the 
mixed 23” series, by W. S. Connor, Contributions to 
probability and statistics (to be published by Stanford 
University Press, 1960). 

“The chain block design, by W. J. Youden and W. S. 
Connor, Biometrics 9, 127 (1953). 

° Experimental designs for duplicate measurements, 
NBS Tech. News Bul. 39, 90 (1955). For further in- 
formation, see Chain block designs with two-way elimi- 
nation of heterogeneity, by J. Mandel, Biometrics 10, 
251 (1954). 

° New experimental designs for paired observations, by 
W. J. Youden and W. S. Connor, J. Research NBS 53, 
191 (1954) ; Partially balanced incomplete block designs 
with two associate classes and two treatments per block, 


by W. H. Clatworthy, J. Research NBS 54, 177 (1955). 


Mathematical Studies of Automatic Mail Sorting 


AN SERIES of mathematical studies for the Post Of- 
fice Department has been conducted by the Bureau 
to help in the design and effective utilization of auto- 
matic mail-sorting equipment. This work comprised 
both a statistical study on the various characteristics of 
mail, including destination, and an analytical study of 
schemes for sorting the mail. Information obtained 
through these studies (currently directed by A. J. Gold- 
man) complements other efforts by the Bureau to speed 
the handling of mail. 

For some time the Post Office Department has been 
seeking mechanized sorting methods to assist its per- 
sonnel in keeping up with the continuing rapid growth 
of letter mail. Because of its specialized experience 
with data processing machines and electronic equip- 
ment, the Bureau has been aiding in various mecha- 
nization problems.‘ For example, last year it arranged 
and supervised a contract for the Rabinow Engineering 
Co. to develop a laboratory prototype machine for 
automatic high-speed sorting. As special techniques 
are required to make optimum use of such equipment, 
the efforts of the Bureau’s mathematicians have also 
been called into play. Speed, efficiency, and economy 
are some of the characteristics which must be optimized 
before automatic sorting equipment can be adopted on 
a wide scale. 


Statistical Study of Characteristics 


A preliminary analysis of the distribution of out- 
going mail was a necessary initial step for estimating 
the feasibility of automatic mail sorting. Special sta- 
tistical techniques were therefore developed for deter- 
mining the flow patterns and physical characteristics of 
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the mail passing through large city post offices. These 
techniques were applied in sampling studies conducted 
in the post offices at San Francisco, Los Angeles, Balti- 
more, and Washington. 

A sampling plan called the “chain ratio” method 
(because its formulas resemble a chain of ratios) was 
successfully used in estimating the distribution of mail 
by destination.” Data of this type had previously been 
gathered by costly procedures involving elaborate 
footage measurements of all the mail passing through 
the sorting system during a single fixed time period, 
such as 24 hr. On the other hand, chain-ratio sam- 
pling, which involved a series of small samples spread 
over a l-week period, required only the counting of 
some 20,000 letters (about 70 ft of mail) daily instead 
of approximately 1,200,000 letters (about 4,000 ft). 

Additional studies * were carried out on the size, 
bulk, and envelope color of mail, on the ratio of hand- 
canceled mail to machine-canceled mail, and on the 
sizes of the top and bottom “clearance spaces” above 
and below the address on an envelope. 

The most important finding of these investigations 
was the uniformity of results from city to city. For 
example, the largest destinations received very nearly 
80 percent of the outgoing mail from San Francisco. 
This was also true of the outgoing mails from Los 
Angeles, and from Baltimore. Again, no significant 
differences with respect to letter-size characteristics of 
canceled mail were found among three cities (San 
Francisco, Los Angeles, and Washington) where 
studies were made. This uniformity indicates that a 
single type of automatic sorting machine can be used in 
large post offices throughout the United States, without 
extensive modifications to allow for local peculiarities. 
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Studies by the applied mathematics 
laboratories have helped in the de- 
sign and effective utilization of auto- 
matic mail sorters. This high-speed 
distributor at the Detroit Post Office 
can handle 54,000 letters an hour 
when all 12 of the operator consoles 
(right) are manned. Start and stop 
controls and speed regulators are 
located inside the raised panel. The 
lights on the end panel indicate the 
sorting rate and signal a warning 
when any of the 279 destination bins 
nears capacity. Photo courtesy Post 
Office Department. 


The methods introduced in this study can be applied 
to other situations in which sampling techniques can 
replace complete enumeration, offering potential sav- 
ings of millions of dollars annually. 


Analytical Study of Sorting Schemes 


The other set of investigations deals with the optimal 
utilization of automatic sorting devices. Several con- 
figurations for this equipment were suggested. One 
(called the Christmas tree scheme) divides, for ex- 
ample, 1,600 destinations into 40 groups of 40. Mail 
is first sorted by group, and then the mail in each group 
is given a second sorting by destination. 

Another proposed configuration (the simple residue 
scheme) involves categorizing each destination as “‘fre- 
quent” or “infrequent” according to the average 
amount of mail addressed to it passing through the 
post office under study. During the first sorting, mail 
to all the infrequent destinations is lumped into a few 
“residue” classifications, and this residue mail then 
receives a second sorting by destination. The point of 
this procedure is that the second sorting brings together 
the residue from several machines, so that the amount 
of mail to even an infrequent destination is increased 
to a point where it can be dealt with more efficiently. 

Analytical methods were developed* for choosing 
among the suggested configurations. The criterion 
was that mail sorting must be carried out at a pre- 
scribed rate for the least cost. When the methods were 
applied to a typical set of data, it was found that the 
Christmas tree scheme had a cost exceeding the mini- 
mum by more than 50 percent. Another proposed con- 
figuration (the straight-forward simplex scheme) led 
to a variable cost which was more than twice the mini- 
mum possible cost. Thus, these configurations are 
unsuitable for large post offices. On the other hand, 
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the simple residue scheme could be brought to within 
8 percent of the minimum cost by optimal selection of 
the breakpoint between the categories of “frequent” 
and “infrequent”. 

A more complicated “hybrid” between two of the 
proposed configurations was developed in which the 
destinations are broken up into three categories of 
“frequent”, “intermediate”, and “infrequent”. By 
this finer subdivision, the cost can be brought within 
4 percent of minimum by optimal choice of the two 
breakpoints separating the categories. 

A separate mathematical study ° concentrated on the 
prevention of overflow in a mail-sorting device devel- 
oped at the Bureau. Earlier investigations elsewhere 
had shown only that a capacity of 120,000 letters (for 
a typical set of data and one of the proposed modes of 
operation) would insure that no overflow occurred. The 
Bureau’s analysis of an identical situation revealed 
that a capacity for only 401 letters would suffice to 
prevent overflow, but that no smaller capacity would 
do so. 


+ Automatic sorting machine for letter mail, NBS Tech. 
News Bul. 42, 156 (1958). 

* A statistician and the post office; a case history in 
operations research, by N. C. Severo, Trans. Conf. of the 
Administrative Applications Division of the Am. Soc. for 
Quality Control 97 (1959); A statistical chain-ratio 
method for determining the distribution of mail by des- 
tination, by NV. C. Severo and A. E. Newman, J. Research 
NBS, 64(C, p. 37 (1960). 

* Some applications of statistical sampling methods to 
outgoing letter mail characteristics, by NV. C. Severo, A. E. 
Newman, S. Young, and M. Zelen, NBS Tech. Note 16. 

* Analytical comparison of suggested configurations for 
automatic mail sorting equipment, by B. K. Bender and 
A. J. Goldman, J. Research NBS 63B, p. 83 (1959). 

° Capacity requirement of a mail sorting device, by 
B. K. Bender and A, J. Goldman, J. Research NBS 62, 
171 (1959); Capacity requirement of a mail sorting de- 
vice, II, by A. J. Goldman, J. Research NBS Sec. 63B, 
p. 79 (1959). 
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Graphical Method for Determining Radio Noise Characteristics 


A? NOISE, particularly atmospheric noise, is the 
limiting factor in radio reception, it is very im- 
portant to be able to predict its characteristics under 
various conditions throughout the radio frequency 
range. Such a prediction, including an estimation of 
the possible interference to a given communications 
system, can be derived from an amplitude-probability 
distribution in which voltage levels are plotted against 
the percentage of time these levels are exceeded. How- 
ever, although the amplitude-probability distribution is 
a particularly useful tool, its direct determination ne- 
cessitates detailed measurements at all frequencies and 
at many locations. The complex equipment and the 
large number of personnel required to carry out these 
observations make continuous routine measurements 
virtually impossible. 

The National Bureau of Standards Boulder Labora- 
tories has found a solution to this problem in an 
empirical graphical method. With this method, it is 
only necessary to measure three statistical moments— 
the average noise power, the average envelope voltage, 
and the average logarithm of the envelope voltage. 
From these quantities, four parameters can be derived 
to describe completely the amplitude-probability dis- 
tribution as plotted on special graph paper designed so 
that a Rayleigh distribution of thermal noise plots as 
a straight line. The curve is the usual plot of voltage in 
decibels against the percentage of time the level is 
exceeded. This approach was developed by W. Q. 
Crichlow, C. J. Roubique, A. D. Spaulding, and W. M. 
Beery of the Bureau’s Central Radio Propagation 
Laboratory. 
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Graph illustrating the effectiveness of the graphical 
method for obtaining the amplitude-probability distri- 
bution of atmospheric noise. The distribution predicted 
from three easily measured moments is seen to agree very 
closely with the distribution measured by complex 
equipment. 
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As the initial step in developing this graphic method, 
it was found that the shape of the curve of over 100 
recorded distributions, when plotted on the special 
paper, could be approximated by two nonparallel 
straight lines connected by the arc of the circle tan- 
gential to the two lines. The lower portion of this 
curve, representing low voltages and high probabilities, 
is composed of many random overlapping events, each 
containing only a small portion of the total energy. 
This portion of the curve approached the Rayleigh 
straight line distribution. The section representing 
very high voltages with low probabilities is, in general, 
composed of nonoverlapping large pulses occurring 
infrequently. This section also was found to be well 
represented by a straight line. 

Four parameters are needed to define this distribu- 
tion: (1) A point through which the Rayleigh line 
passes; (2) a point and (3) a slope for the high- 
voltage, low-probability line; and (4) a parameter of 
some kind to determine the radius of the circular arc. 
Since four parameters are needed and only three 
statistical moments are measured, a dependency be- 
tween two of the parameters was established. From 
experimentally measured distributions at eight fre- 
quencies and various bandwidths, a linear relationship 
was found between the radius of the circle and the slope 
of the high-voltage, low-probability line. 

The experimentally measured distributions were 
then numerically integrated to obtain the three 
measurable moments, so that the relationship between 
the three moments and the four parameters could be 
determined. It is this relationship which allows 
measured moments to be used to plot the amplitude- 
probability distribution. 

As a check on the accuracy of this method, the 
numerically integrated moments were compared to 
moments which had been recorded at the same time as 
the distributions. Approximately 60 sets of distribu- 
tion measurements were analyzed to determine the 
error, which was largely introduced by the recorder. 
The standard deviation in the rms voltage values was 
about 1.1 db, with an average error of —0.4 db. 
Ninety percent of the values were within 2 db of the 
mean. 

This method of obtaining an amplitude-probability 
distribution from the three measured moments results 
in a distribution which is valid only for the bandwidths 
in which the moments were measured. A study is now 
in progress to develop an accurate method of band- 
width conversion. 


1 Determination of the amplitude-probability distribu- 
tion of atmospheric radio noise from statistical moments, 
by W.Q. Crichlow, C. J. Roubique, A. D. Spaulding, and 
W.M. Beery, J. Research NBS 64D p. 49 (1960). 
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Mathematical Research on Reliability Prediction 


4 pee IMPROVEMENT of reliability is one of the 
major problems facing current technology. Elec- 
tronic computers, automatic control systems, and 
satellite instrumentation must operate with satisfactory 
reliability over relatively long periods of continuous 
or intermittent use. On the other hand, the electronic 
systems in guided missiles are required to withstand a 
long “shelf-life” yet perform with high reliability dur- 
ing their very brief period of use. 

Because of the complexity of such systems, many 
unknown and uncontrollable factors influence their 
performance. Thus, some aspects of system perform- 
ance have come to be treated as probabilistic phe- 
nomena. For quantitative evaluation, the reliability 
of a system may be defined as the probability of satis- 
factory performance under clearly specified conditions, 
including time of operation. The evaluation and pre- 
diction of system reliability accordingly require the 
development and use of probabilistic and statistical 
tools. 

The Bureau is studying problems in mathematical 
statistics and applied probability theory that arise in 
evaluating the reliability of complex systems. Empha- 
sis is placed on alternative approaches to the construc- 
tion of mathematical models for the organization of 
information on system design and component perform- 
ance. This work, by J. R. Rosenblatt and M. Zelen of 
the Statistical Engineering Laboratory, is comple- 
mented by the development of statistical techniques for 
analyzing experimental data on the reliability of com- 
ponents. The long-range goal of this program is to 
make available a spectrum of approaches to model 
construction and a variety of statistical techniques, 
with accompanying criteria for selecting the most 
suitable method for a particular application. 


System Reliability 


Evaluation or prediction of system reliability, espe- 
cially at the research and development stage, depends 
on systematic utilization of engineering design infor- 
mation, data on component performance, and data on 
the conditions under which the system will be used. 
Studies at the Bureau have led to the formulation of 
a general framework for mathematical representation 
of causal or statistical interdependence of subsystems 
in a complex system.1. This work provides a basis for 
the detailed treatment of the relation of system per- 
formance to properties of system design and component 
performance. 

Within the general framework, various possible 
forms of relationships among subsystems have been 
studied. It was shown, for example, how analysis of 
system design and operating environments may lead to 
separation of environmental effects from other factors. 
When this can be done, subsystems may be found to 
be (statistically) conditionally independent under 
fixed environmental conditions. When conditional 
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independence holds, the calculation of a reliability pre- 
diction will be simplified. 

The classification of modes of failure is one of the 
preparatory steps in constructing a mathematical 
model for reliability prediction. The definition of 
different phases of the storage and operating life of the 
equipment under consideration is another step. The 
Bureau has studied the way in which reliability predic- 
tions may be affected by the choice of classifications 
and definitions.” These choices are sometimes forced 
by the form of available data, so that an existing data 
collection program may to a large extent determine the 
relevance and meaningfuless of a reliability predic- 
tion. On the other hand, explicit examination of 
alternative possibilities for classification and definition 
may help in designing a program for obtaining more 
useful subsystem and component data. 


Component Reliability 


In order to assess the reliability of a system, it is 
necessary to evaluate the reliabilities or lifetimes of 
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The Bureau computed these operating characteristic 
(OC) curves for a sequential life test procedure based 
on exponential distribution theory. If the distribution 
of failure times is not in fact exponential, the OC curves 
may be quite different from the curve shown for the 
exponential case—as is the OC curve obtained when the 
failure times follow a certain Weibull distribution. In 
this particular case, component lots having an average 
lifetime of 1,000 hr would be accepted with probability 
0.90 if the life distribution were exponential, whereas if 
the life distribution were of the Weibull form, the ac- 
ceptance probability would be 0.90 for lots having an 
average life as small as 760 hr. Similarly, lots haying 
an average life of 500 hr would be accepted with prob- 
ability 0.10 if the underlying distribution were expo- 
nential but would be accepted with probability 0.00 if the 
underlying distribution were of the Weibull form. This 
work is part of a Bureau study on the “robustness” of 
statistical life-testing procedures for acceptance sampling. 
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the various components. Data for this purpose are 
provided by “life tests’—experiments in which the 
lifetime of a particular type of part is measured. 

The Bureau has developed techniques for analyzing 
life test data obtained under several sets of combined 
environmental conditions. The techniques include 
statistical tests of significance which evaluate the im- 
portance of environmental factors, singly or in com- 
bination, and procedures which allow confidence 
intervals to be calculated for the relative effect of each 
experimental factor on mean life. A statistical proce- 
dure was derived which gives some indication of the 
meaningfulness of “accelerated” life tests, that is, life 
tests conducted under extreme environmental 
conditions. 

The applicability of statistical procedures for life 
tests, when the underlying assumptions are not satisfied, 
has also been investigated.* All statistical techniques 
are based on assumptions concerning the underlying 
probability distributions of the phenomena observed. 
Techniques which can be used when their special as- 
sumptions are not exactly satisfied are called “robust”. 

A number of widely used statistical life testing proce- 


dures are based on the assumption that the distribution 
of lifetime has the exponential form. The Bureau found 
that some of these procedures have very poor robust- 
ness. Further work on the robustness of life testing 
techniques is planned; it will include both theoretical 
studies and computer simulation of life tests. 

Another type of component experiment also presents 
problems in theoretical statistics. This is the “aging 
experiment”, in which life histories of component per- 
formances, as well as length-of-life data, are collected. 
Work now in progress deals with statistical techniques 
that may be used for the analysis of aging experiments. 


*On prediction of system performance from informa- 
tion on component performance, by J. R. Rosenblatt, 
Proc. Western Joint Computer Conf. 85 (February 1957). 

?On prediction of system behavior, by J. R. Rosen- 
blatt, Proc. NYU-Industry Conference on Reliability 
Theory 39 (June 1958). 

* Factorial experiments in life testing, by M. Zelen, 
Technometrics 1, 269 (1959). 

* Are life testing procedures robust? by M. Zelen and 
M. C. Dannemiller, Proc. Sixth National Symposium on 
Reliability and Quality Control in Electronics (to be 
published) (January 1960). 


Mathematical Studies in Plasma Physics 


HE BUREAU, under the sponsorship of the Air 

Research and Development Command, has initiated 
a theoretical research study aimed at solving basic 
problems in plasma physics. This work, by C. M. 
Tchen of the mathematical physics group, forms a part 
of a comprehensive Bureau-wide program to provide 
basic data and theory for plasma physics and astro- 
physics. As plasma is a unique state of matter with 
properties depending strongly on its dynamical con- 
figuration, the present project will emphasize dynam- 
ical aspects. That is, it will concentrate on the flow 
of an ionized gas through magnetic and electric fields. 

In the past few years, the field of theoretical plasma 
dynamics and magnetohydrodynamics has grown 
rapidly, creating a wide variety of research problems. 
Moreover, the recent developments of space technology 
and cosmic flight have focused attention on the dy- 
namic behavior of cosmic gases, their interaction with 
magnetic fields, and the motion of bodies—stars, plan- 
ets, and artificial satellites—in a cosmic medium. Be- 
cause of the increasing importance of the field and the 
present lack of basic data, the present research program 
was undertaken. 

A large part of this theoretical work will be devoted 
to establishing sound foundations for plasma theory. 
The Bureau is presently developing a quantum theoreti- 
cal generalization using the methods of quantum colli- 
sion theory. Relativity considerations will later be 
included to arrive at a general formulation, which will 
be applicable to cases where electrons attain relativistic 
energies. 

Another study being planned will emphasize the 
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transport properties of plasmas and the statistically 
determined collective behavior on which such prop- 
erties depend. An important phase of this work deals 
with the effect of the electric field on the motion of the 
charged particles which make up the plasma. A sys- 
tematic study of the correlations in a plasma seems to 
be a natural and promising starting point for develop- 
ing a theory of the effective electric field. 

Studies to describe the energy dissipated through 
radiation during a rise in plasma temperature are also 
being pursued. An elementary theory of brems- 
strahlung loss which can be applied to plasma is pres- 
ently available. However, its limits of applicability 
still have to be investigated from the point of view of 
a more general theory—one which will take into ac- 
count relativistic considerations. A theory of the 
radiation emitted by charged particles revolving in a 
magnetic field is in principle available both for non- 
relativistic and relativistic energies. The Bureau in- 
tends to investigate its applicability to plasma con- 
ditions. 

Another goal of the program is the extension of the 
magnetohydrodynamic equations to take into account 
a nonscalar tensor pressure. ‘This will allow the 
hydrodynamic equations, modified to include the 
interaction between the motion of the fluid and the 
magnetic field, to be applied to gas discharges and 
certain astrophysical phenomena. 

With these and other theoretical studies, the Bureau 
hopes to provide the foundations for further develop- 
ments in fusion research, the physics of space, and the 
study of stellar and planetary atmospheres. 
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Research on Satellite Orbis 


N WORK done for the Air Force, the Bureau has 
devised a simple, direct approach for including the 
effects of the earth’s oblateness, or flattening at the poles, 
in the calculation of satellite orbits.t Based on an anal- 
ysis by J. P. Vinti of the mathematical physics group, 
this approach applies in cases where drag is negligible. 
Although results are still being worked out, it promises 
to eliminate many of the difficulties encountered in 
applying perturbation theory to the problem of calcu- 
lating the effects of oblateness. Furthermore, it can 
be used even when the orbital inclination approaches 
the critical value of 63.4 deg, where special considera- 
tions might otherwise be necessary. 

A single passage of an artificial satellite about the 
earth is approximately an ellipse. If the earth had 
true spherical symmetry, then a true ellipse, relative 
to an inertial frame of reference, would be described 
under certain conditions. These conditions are that 
the perigee * be high enough so that the drag be negli- 
gible and that the apogee * be low enough so that the 
moon’s action also be negligible. However, because 
of the oblateness of the earth, the ascending node of 
the orbit has a westward secular motion of 3 or 4 deg 
a day, along with a small periodic motion. Oblateness 
also gives rise to a motion of the line of apsides (the 
major axis of the “ellipse’”) relative to the line of nodes 
(the line from the center of the earth to the ascending 
node). 

Previous investigators have applied perturbation the- 
ory to find the details of these effects of oblateness, 
beginning either with an unperturbed orbit corre- 
sponding to true spherical symmetry or with an inter- 
mediate orbit that allows for part of the force 
produced by oblateness. These intermediate orbits can 
be found exactly in terms of quadratures, but they al- 
low only partly for oblateness and correspond to gravi- 
tational potentials that do not satisfy Laplace’s 
equation. 

The Bureau’s approach utilizes oblate spheroidal 
coordinates to take advantage of the geometrical sym- 


Some of the effects of the earth’s 
oblateness, or flattening at the poles. 
Because the earth is not a perfect 
sphere, an artificial satellite does not 
describe a true ellipse. The ascend- 
ing node of the orbit has a westward 
secular motion (left). In addition, 
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oblateness gives rise to a motion of PLANE OF 
the line of apsides (center). The EQUATOR 
combined effect is illustrated at 

right. Previous treatments have LINE OF 
used perturbation theory in calcu- NODES 


lating these and other effects caused 
by oblateness. The Bureau recently 
devised a more direct approach to 
the problem which takes into ac- 
count the full oblateness of the 
earth at the start of the calculation. 


Photo courtesy of NASA, 
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metry of the problem. This approach involves finding 
the most general solution of Laplace’s equation, non- 
singular outside the earth, that results in separability 
of the Hamilton-Jacobi equation. The resulting grav- 
itational potential contains certain parameters which 
can be adjusted to fit the 1/r term and the second 
harmonic in the empirically known expansion in spher- 
ical harmonics. The higher even harmonics are then 
fixed, but it turns out that the fourth harmonic thus 
obtained agrees, within the limits of observational er- 
ror, with the recent value found by D. G. King-Hele 
and and R. H. Merson * by analysis of data from ob- 
served satellite orbits. The corresponding “inter- 
mediate orbit”, which can also be found exactly in 
terms of quadratures, allows so accurately for oblate- 
ness that it may be considered to be close to a final 
orbit, so far as even harmonics are concerned. 

For investigating the effects of other factors, such as 
equatorial asymmetry represented by odd harmonics, 
local gravitational anomalies represented by tesseral 
harmonics, and drag, the solution can then be used as 
the starting point for perturbation theories. These 
other effects are all expected to be small compared with 
that of oblateness for orbits that are essentially above 
the atmosphere. 


* New method of solution for unretarded satellite orbits, 
by J. P. Vinti, J. Research NBS 63B (in press). 

*The point in the orbit of a satellite nearest to the 
earth’s center. 

* The point in the orbit of a satellite farthest from the 
earth’s center. 

*A new value for the earth’s flattening devised from 
measurements of satellite orbits, by D. G. King-Hele and 
R. H. Herson, Nature 183, 881 (1959) ; Recent studies 
of satellite orbits, by E. C. Cornford, D. G. King-Hele, 
and R. H. Merson (Institute of the Aeronautical Sciences 
Paper No. 59-141). Note added in press: A later value 
for the fourth harmonic, by O’Keefe, Eckels, and Squires, 
Astronomical Journal 64, 245 (1959), does not appear 
to be compatible with the value here mentioned. Such 
a change, if confirmed, could produce appreciable effects 
only near the critical inclination. 
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Complete sequences and approximations in normed linear 
spaces. P. Davis and K. Fan. Duke Math. J. 24, 183 

(1957). 

Numerical experiments in potential theory using orthonormal 
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46, 12 (1956). 
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The mean deviation of the Poisson distribution. E. L. Crow. 
Biometrika 45, 556 (1958). 
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